** Figure 5

clear all
set more off 
capture log close
 
* main path to data
global path "~/Dropbox/Paraguay/1_data/7_replication_dataverse"
capture cd "$path"	

*figures= analysis figures
global figures "$path/figures"

* main panel
use "$path/data/main.dta", clear

xtset id year

// Generate indicator variables for years 1954 to 2007
forval year = 1954/2007 {
    generate y`year' = (year == `year')
}

// Generate the alt_1989 variable and set it to 0
generate alt_1989 = 0


** Figure 5a

set scheme s1mono

set more off
reghdfe tmh_ha_log ///
	y1954#c.zeigen y1955#c.zeigen y1956#c.zeigen y1957#c.zeigen y1958#c.zeigen ///
	y1959#c.zeigen y1960#c.zeigen y1961#c.zeigen y1962#c.zeigen y1963#c.zeigen ///
	y1964#c.zeigen y1965#c.zeigen y1966#c.zeigen y1967#c.zeigen y1968#c.zeigen ///
	y1969#c.zeigen y1970#c.zeigen y1971#c.zeigen y1972#c.zeigen y1973#c.zeigen ///
    y1974#c.zeigen y1975#c.zeigen y1976#c.zeigen y1977#c.zeigen y1978#c.zeigen ///
    y1979#c.zeigen y1980#c.zeigen y1981#c.zeigen y1982#c.zeigen y1983#c.zeigen ///
    y1984#c.zeigen y1985#c.zeigen y1986#c.zeigen y1987#c.zeigen y1988#c.zeigen ///
    alt_1989 ///
    y1990#c.zeigen y1991#c.zeigen y1992#c.zeigen y1993#c.zeigen y1994#c.zeigen ///
    y1995#c.zeigen y1996#c.zeigen y1997#c.zeigen y1998#c.zeigen y1999#c.zeigen ///
    y2000#c.zeigen y2001#c.zeigen y2002#c.zeigen y2003#c.zeigen y2004#c.zeigen ///
    y2005#c.zeigen y2006#c.zeigen y2007#c.zeigen ///
    zeigen [pweight=p_weights] if sample_weight == 1, absorb(year id) cluster(id) nocons

coefplot, vertical omitted drop(zeigen y1954-y1988 y1990-y2007 ///
    1.y1954#c.zeigen 1.y1955#c.zeigen 1.y1956#c.zeigen 1.y1957#c.zeigen 1.y1958#c.zeigen ///
    1.y1959#c.zeigen 1.y1960#c.zeigen 1.y1961#c.zeigen 1.y1962#c.zeigen 1.y1963#c.zeigen ///
    1.y1964#c.zeigen 1.y1965#c.zeigen 1.y1966#c.zeigen 1.y1967#c.zeigen 1.y1968#c.zeigen ///
    1.y1969#c.zeigen 1.y1970#c.zeigen) ///
    yline(0, lcolor(red)) ///
    xlabel(1 "1971" 4 "1974" 7 "1977" 10 "1980" 13 "1983" 16 "1986" 19 "1989" 22 "1992" ///
    25 "1995" 28 "1998" 31 "2001" 34 "2004" 37 "2007", angle(45)) ///
    xtitle("Year", size(3.7) margin(small)) ///
    ytitle("Effect relative to 1989", size(3.7) margin(small)) ///
    ciopts(recast(rcap)) ///
    yscale(range(-0.05 0.2)) ///
    ylabel(-0.05(0.05)0.2, angle(0)) ///
    xline(19, lpattern(dash))
	
	
	

graph export "$figures/figure5a.pdf", replace

	

** Figure 5b


set more off

reghdfe tmh_e ///
    y1954#c.zeigen y1955#c.zeigen y1956#c.zeigen y1957#c.zeigen y1958#c.zeigen ///
	y1959#c.zeigen y1960#c.zeigen y1961#c.zeigen y1962#c.zeigen y1963#c.zeigen ///
	y1964#c.zeigen y1965#c.zeigen y1966#c.zeigen y1967#c.zeigen y1968#c.zeigen ///
	y1969#c.zeigen y1970#c.zeigen y1971#c.zeigen y1972#c.zeigen y1973#c.zeigen ///
    y1974#c.zeigen y1975#c.zeigen y1976#c.zeigen y1977#c.zeigen y1978#c.zeigen ///
    y1979#c.zeigen y1980#c.zeigen y1981#c.zeigen y1982#c.zeigen y1983#c.zeigen ///
    y1984#c.zeigen y1985#c.zeigen y1986#c.zeigen y1987#c.zeigen y1988#c.zeigen ///
    alt_1989 ///
    y1990#c.zeigen y1991#c.zeigen y1992#c.zeigen y1993#c.zeigen y1994#c.zeigen ///
    y1995#c.zeigen y1996#c.zeigen y1997#c.zeigen y1998#c.zeigen y1999#c.zeigen ///
    y2000#c.zeigen y2001#c.zeigen y2002#c.zeigen y2003#c.zeigen y2004#c.zeigen ///
    y2005#c.zeigen y2006#c.zeigen y2007#c.zeigen ///
    zeigen [pweight=p_weights] if sample_weight == 1, absorb(year id) cluster(id) nocons

	
	
coefplot, vertical omitted drop(zeigen y1954-y1988 y1990-y2007 ///
    1.y1954#c.zeigen 1.y1955#c.zeigen 1.y1956#c.zeigen 1.y1957#c.zeigen 1.y1958#c.zeigen ///
    1.y1959#c.zeigen 1.y1960#c.zeigen 1.y1961#c.zeigen 1.y1962#c.zeigen 1.y1963#c.zeigen ///
    1.y1964#c.zeigen 1.y1965#c.zeigen 1.y1966#c.zeigen 1.y1967#c.zeigen 1.y1968#c.zeigen ///
    1.y1969#c.zeigen 1.y1970#c.zeigen) ///
	    yline(0, lcolor(red)) ///
    xlabel(1 "1971" 4 "1974" 7 "1977" 10 "1980" 13 "1983" 16 "1986" 19 "1989" 22 "1992" ///
    25 "1995" 28 "1998" 31 "2001" 34 "2004" 37 "2007", angle(45)) ///
    xtitle("Year", size(3.7) margin(small)) ///
    ytitle("Effect relative to 1989", size(3.7) margin(small)) ///
    ciopts(recast(rcap)) ///
    yscale(range(-0.01 0.03)) ///
    ylabel(-0.01(0.01)0.03, angle(0) format(%9.2f) labsize(small)) ///
    xline(19, lpattern(dash))

graph export "$figures/figure5b.pdf", replace
	
	
	
